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ABSTRACT 

Herschel PACS and SPIRE images have been obtained of NGC 6720 (the Ring Nebula). This is an evolved planetary nebula with a central star 
that is currently on the cooling track, due to which the outer parts of the nebula are recombining. From the PACS and SPIRE images we conclude 
that there is a striking resemblance between the dust distribution and the H2 emission, which appears to be observational evidence that H2 forms 
on grain surfaces. We have developed a photoionization model of the nebula with the Cloudy code which we used to determine the physical 
conditions of the dust and investigate possible formation scenarios for the H2. We conclude that the most plausible scenario is that the H2 resides in 
high density knots which were formed after the recombination of the gas started when the central star entered the cooling track. Hydrodynamical 
instabilities due to the unusually low temperature of the recombining gas are proposed as a mechanism for forming the knots. H2 formation in the 
knots is expected to be substantial after the central star underwent a strong drop in luminosity about one to two thousand years ago, and may still 
be ongoing at this moment, depending on the density of the knots and the properties of the grains in the knots. 
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1. Introduction 

Grains play an important role in many environments, includ- 
ing planetary nebulae (PNe), because of extinction, photoelectric 
heating, their influence on the charge and ionization balance of 
the gas, as catalysts for grain-surface chemical reactions (e.g. H2 
formation), and as seeds for freeze-out of molecules. Previous 
satellite missions such as IRAS, ISO, Spitzer, and AKARI have 
allowed us to study the dust in PNe, but unfortunately the an- 
gular resolution of these instruments was too low to get detailed 
information on the spatial distribution of the dust. This has now 
changed with the launch of Herschel, which allows us to study 
the spatial structures in unprecedented detail. In this paper we 
will do this for NGC 6720 (M57, the Ring nebula) to study H2 
formation. NGC 6720 is an evolved, oxygen-rich bipolar nebula 
seen nearly pole-on. The nebula is optically thick to ionizing ra- 
diation in most dire ctions, but optically thin in the polar regions 
JO'Dell et al.ll2007l her eafter OD07). Molecules such as H2 and 
CO h ave been detected jBeckwith et aI.lll978l: lHuggins & Healvl 
Il986l) . The central star has exhausted hydrogen shell burning and 
is now on the cooling track. As a result the outer halo is recom- 
bining and re-ionization of the innermost recombined material 
due to expansion of the nebula has just started (OD07). This ob- 
ject is very similar to the Helix nebula, which seems to be further 
advanced along the same evolutionary path (OD07). In Sect. |2] 
we will describe the observations, and in Sect. [3] we will discuss 
various scenarios for the formation of H2. 



* Herschel is an ES A space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA. 
** email: p.vanhoof@oma.be 



2. Observations 

The images of NGC 6720 presented in this paper were obtained 
with the PACS and SPIRE instruments on board the Herschel 
satellite (Pilbratt et al. 2010) on 2009- 10- 10 and 2009- 10-06, re- 
spectively, as part of the Science Demonstration Phase of the 
Mass-loss of Evolved StarS (MESS) guaranteed time key pro- 
gram (Gi'oeneweg en et al., in preparati on). The PACS instrument 
is described by Poglits ch et al.l (^010). The SPIRE instrument, 
its in-orbi t performance, and its scientific capabilities are de- 
scribed bv lGriffin et al.l ( 1201 Oh . and the SPIRE a stronomical cal- 
ibration methods and accuracy are outlined by ISwinvard et al.l 
(2010). The reduction of our data is described in Appendix lAl 
We also determined the total flux in each of the images. This is 
described in Appendix |B] Below we will discuss these images 
in detail and compare them to the H2 morphology of the nebula. 
To this end we used a ground-based H2 2. 12 jum image, obtained 
with the Omega2000 camera on the 3.5-m Zeiss telescope at the 
Calar Alto Observatory. 

In Fig. [T] we present the PACS and SPIRE images and com- 
pare them to the Calar Alto H2 image. The morphology of the 
dust and the H2 emission shows a striking resemblance even 
in small details. This appears to be observational evidence that 
H2 forms on grain surfaces in an astrophysical environment. 
Such eviden c e is v ery rare. A similar result was obtained by 
iHabart et al.l (l2003h for the p Ophiuchi molecular cloud, sug- 
gesting that H2 forms on PAH surfaces. The presence of PAHs 
in NGC 6720 cannot be fully excluded, but seems very unlikely. 
Accordingly our observations are the first indication for H2 for- 
mation on oxygen-rich dust grains in an astrophysical environ- 
ment to our knowledge. 
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Fig. 1. NGC 6720 in five photometric bands. Top row from left to right: H2 2. 12 fim, PACS 70 fim, PACS 160 yum. Bottom row from 
left to right: SPIRE 250 fim, SPIRE 350 yum, and an overlay of the Calar Alto H2 contours on the PACS 70 fim image. The H2 image 
is not flux calibrated. The maps have standard orientation (N to the top, E to the left). 



We can divide the H2 image into 3 regions: the inner ring, 
with a semimajor axis of ~ 45", an inner set of arcs with a radius 
of ~ 70" (the inner halo), and an outer set of arcs with a radius 
of ~ 1 10" (the outer halo). There is also fainter emission outside 
of these outer arcs, but this is not very easy to see on the printed 
map. 

Comparing the H2 and the sub-mm maps it is clear that the 
size and shape of the inner ring is reproduced on all images. The 
same can be said for the inner halo, which can be seen clearly 
on all maps. In the NW part of the inner halo, the H2 emission 
region is broader than in the SE — this too is the case on the 
SPIRE images, and less clearly so on the PACS images. Finally, 
the outer halo on the H2 image is also visible on the sub-mm 
images — as extended emission on the SPIRE images and as a 
faint circular ring shape on the PACS images (on the 70 //m im- 
age most clearly). On the zoomed PACS 70yum with overlaid H2 
contours (Fig.[T] lower right panel) the clear correspondence be- 
tween the optical and sub-mm emission in the ring and the inner 
halo is highlighted. 



3. The origin of the H2 

From the ground-based H2 images it is clear that most of the H 2 
resides in high density knots in the inner ring dSpeck et al.l2003h . 
In this section we will investigate the origin of this H2. The halo 
also shows H2 emission, which has a very different morphol- 
ogy though. The latter will not be discussed here. It is clear 
that H2 was formed in the dense AGB wind. In the post-AGB 
phase three different scenarios will be investigated: 1) this H2 
survived in the ionized region, 2) this H2 survived in the knots. 



Table 1. Parameters of the Cloudy model of NGC 6720. 



Teff (kK) 


134.7 


e(C) 


8.20 


U (Lo) 


239. 


6(N) 


7.89 


Hn (mpc) 


56. 


e(0) 


8.50 


r„„, (mpc) 


175. 


6(Ne) 


7.93 


log(«H) (cm-3) 


2.60 


e(S) 


6.42 


T, (kK) 


11.97 


f(Cl) 


4.94 


log(ne) (cm-3) 


2.62 


6(Ar) 


6.19 


log(r) 


-2.21 


6(Fe) 


5.18 


e(He) 


11.00 


D(pc) 


740. 



which formed before the gas was ionized and 3) this H2 was de- 
stroyed and then was formed again later inside the knots when 
they formed. These scenarios were already investigated for the 
Helix nebula bv lMatsuura et al.l (120091) . In order to determine the 
physical conditions in the nebula we created a photoionization 
model using a prerelease of version CI 0.00 of the photo ioniza- 
tion co de Cloudy (revision 3862), last described by Ferland et alj 
(Il998h . The details of the modeling are described in AppendixICl 
the resulting model parameters are summarized in Table [T] 



3.1. Can H2 survive in tine ionized region? 

For the Helix nebula. Matsuura et al.l (l2009l) conclude that "part 
of the H2 is primordial, i.e. formed during the AGB phase, and 
survived the ionization of the nebula". This conclusion is based 
on the work of Aleman & Gruenwald (2004, hereafter AG04). 
We recreated the models of AG04 with Cloudy. This code is ca- 
pable of creating an equilibrium model of the ionized region as 
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well as the PDR and the molecular r egions. It includes st ate- 
of-the-art code for modeling grains ("van Hoof et alJ [20041) as 
well as th e chemistry network (Abeletal. 2005) and the H2 
molecule dShaw et alJ i2005). We created spherically symmet- 
ric, one-dimensional models with the parameters of the standard 
model of AG04. Our H2 model includes all ro-vibrational levels 
of the ground electronic state, as well as the six electronically 
excited states that are coupled to the ground state by permit- 
ted electronic transitions. Self-shielding of the H2 molecule was 
calculated self-consistently using detailed radiative transfer (in- 
cluding line overlap) of each of the roughly half a million lines 
of the molecule. 

First we needed to define the outer radius of the ionized re- 
gion. In PNe that is not straightforward. Due to the v"^ depen- 
dence of photoionization cross section, the high energy tail of the 
stellar spectrum will propagate through the ionized region into 
the PDR and cause hydrogen to stay partially ionized there. Such 
PDRs are sometimes referred to as XDRs. The transition from 
fully ionized gas to neutral gas will be gr adual and the transition 
layer w ill have a considerable thickness. lAleman & Gruenwaldl 
(|2004|) stopped their calculations when n(H^)/nH - 10""* was 
reachecQ. Cloudy by default stops when the electron tempera- 
ture drops below 4000 K, based on the notion that at lower tem- 
peratures the forbidden optical emission lines will no longer be 
excited effectively and the stopping radius will coincide with the 
outer radius in optical images of the nebula. The choice of the 
stopping criterion is crucial, but is unfortunately not discussed 
by AG04. The H2 fraction Rm (Eq. 7 in AG04) continually rises 
as the value of n{H^)/nu where the calculation is stopped, is 
lower. This is illustrated in Fig.|2] This clearly shows that there 
is no obvious choice for the stopping criterion: the deeper one 
integrates, the higher the value of R^ becomes. This is hardly 
surprising: deeper layers are better shielded from UV radiation 
and the H2 abundance will be higher. The local maximum of the 
H2 abundance that AG04 report around n(H^)/nu = 0. 1 is barely 
reproduced in the Cloudy model and has virtually no effect on 
the overall value of Rm (see Fig. [3]). Most of the H2 present 
in this model will be in regions near the outer radius. The H2 
weighted average ionization of hydrogen is slightly more than 
10"^. So a more precise statement would be that a small fraction 
of H2 can survive (or better: is constantly destroyed and formed 
again) in the transition region between the ionized region and 
the molecular region. The expectation of AG04 that using a de- 
tailed treatment of H2 self-shielding could increase the value of 
Rfa by up to 2 dex could not be confirmed with Cloudy, at least 
not for the standard model of AG04. Our value of Rm is actu- 
ally slightly lower than reported by AG04. The standard model 
of AG04 is roughly appropriate for a low-mass central star on 
the horizontal track. To model a high-mass central star, we also 
created a "standard4" model which is identical to the standard 
model except that hh = lO'* cm"^ and L, - 10^ Lq. The results 
are shown in Figs. |4]and|5] In this model, the value for Rm is 
somewhat lower than in the standard model. So the conclusion 
is that a small amount H2 can exist in the ionized region, mostly 
in the transition zone towards the PDR, and that this amount is 
not sufficient to explain the observed H2 emission in the inner 
ring. 



3.2. Can H2 survive in knots? 

When the photoionization of the AGB shell starts at the moment 
when a planetary nebula is born, the molecules in the circumstel- 
lar shell will be very quickly destroyed unless they can somehow 
be shielded from the ionizing radiation. Dense knots that formed 
during the AGB phase can provide this environment and seem 
a very attractive explanation for the existence of H2 in evolved 
PNe such as NGC 6720. It is clear that dense knots can shield 
the material in the center from FUV flux and allow H2 molecules 
to s urvive. But it is also obvious from the images presented in 
e.g. iMatsuura et all ( l2009i) that the FUV flux will interact with 
the outer layers of the knot. This interaction causes material to 
be advected off the knot. During that process the material will 
be heated, causing H2 to shine. This provides a na tural explana- 
tion f or the H2 spectrum seen in the Helix nebula (iHennev et alJ 
2007). This process will gradually erode the knots, and eventu- 
ally they will be destroyed. We will now make an estimate of the 
lifetime of a knot in the ionized region. 

We adopt a value of 1.5x10'^ M p for the mass of the knots 
(iMeixner et al.ll2005l; IMatsuura et"aLl l2009). We furthermore as- 
sume that the knots have a typical diameter of 0.5". Combined 
with a distance of 740 pc ( O'Dell et al. 20 09) this gives a phys- 
ical diameter of 5.54x1 0^° km (37 AU), in agreement with the 
typical size stated in .O'Dell et all ( l2003h . Assuming that 70% 
of the mass is hydrogen, this yields an average number density 
of hydrogen of nn =1.4x10^ cm ■'. The density in the center of 
th e knot should be highe r, in agreement with the m easurements 
bv lMeaburnetal.1 (Il998h and lHuggins et al.l (l2002h . We assume 
spherical astronomical silicate grains to be present in the knot, 
with a standard ISM size distribution (Mathis et al. 1977) and a 
dust-to-gas mass ratio equal to the standard model of NGC 6720. 

iHennev et aH (l2007l) state that their models C06 and A06 
matched the observations of the Helix nebula best for the outer 
and inner knots, respectively. Taking the parameters for the ad- 
vection flow presented in Table 1 of Hennev et al. (2007), one 
can calculate that the erosion rate would be approximately 10"'" 
and 10"^ Moyr ' for the C06 and A06 model, respectively. 
Combining this with the adopted mass of the knot yields a sur- 
vival time of 15,000 to 150,000 yr. However, one should realize 
that during most of PN evolution the central star luminosity was 
much higher than assumed here (by upto 2 dex, the central star 
of the Helix nebula has a luminosity of 120 L©), and also that 
the nebula was much more compact during the early stages of 
evolution causing the knots to be much closer to the central star 
The will cause the erosion to be much faster early on and hence 
the survival time to be much shorter The erosion rate will scale 
at least with the square-root of the EUV flux when the advec- 
tion flow is in the recombination-dominated regime. So a con- 
servative estimate is that the survival time i s at least a fa c tor 10 
shorter than stated above, and likely more. lO'Dell et al.l (l2007l 
120091) state that the kinematic age of the Ring nebula is 7000 yr, 
making the survival of the knots inside the ionized region prob- 
lematic. Note that this argument also applies to models where 
the knots formed due to instabilities at the ionization front dur- 
ing the onset of ionization. More detailed modeling of the knots 
during the high-luminosity phase is warranted though to reach a 
more definitive conclusion. 



' The number 10"'' is stated in tfie text, but an inspection of their 
Fig. 1 suggests that they actually used 10"^^ . In this section we will adopt 
10""* as the stopping criterion. 



3.3. Can H2 be formed again after ionization? 



lO'Dell et al.l (l2007h argue that knot formation only starts after 
the central star has entered the cooling track and the nebular ma- 
terial starts recombining. This argument is based on the fact that 
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no similar features are seen in HST images of younger PNe. We 
propose a new explanation in which the recombining gas cools 
very quickly, much fast er than the recombination proceeds. This 
is indicated by Fig. 4 in lFerland & Truranl(ll98lh . but these cal- 
culations clearly need to be redone for conditions that are more 
appropriate for NGC 6720. The fast cooling would result in gas 
that is still ionized (and thus produces abundant recombination 
radiation) but has very little thermal pressure support. The radi- 
ation pressure of the recombination radiation on the dust and/or 
gas may then cause the medium to become unstable and frag- 
me nt into many g lobule s . 

lO'Dell et al.l (120071) also argued that the main gas phase of 
NGC 6720 was fully ionized during the high-luminosity phase, 
so that any H2 present in the gas would have been destroyed very 
quickly (on timescales less than a year). Consequently any H2 
present in the knots must have formed after the formation of the 
knots itself started and they became dense enough to shield the 
molecules from the UV light. To investigate this scenario further, 
we created Cloudy models of the knots using revision 3862 of 
the experimental "newmole" version of the code in order to cope 
with the very steep density gradient at the outer edge of the knot. 
As was the case with the models described in Sect. 13. II we use a 
full and self-consistent H2 model. The models are based on our 
standard model of NGC 6720, with a modified density law that 
mimics a dense knot just outside the ionized region. The density 
is constant until the electron temperature drops below 4000 K, 
then the density will rise steeply towards the core density of the 
knot with an assumed scale-length of 2x10'^ km. 

We computed eight Cloudy models with log nu - 4 (0.5) 7.5 
to cover the full range of plausible densities in the center of 
the knot. From the models we extracted the formation (solid 
lines) and destruction timescales (dotted lines) of H2. The re- 
sults are shown in Fig. |6] for selected models. The Cloudy mod- 
els show that a minimum core density of roughly 10^ cm"-' is 
needed to get sufficient shielding from the U V radiation for sub- 
stantial H2 formation in the knot. lGould & Salpeter ( 196 3) have 
shown that H2 formation is dominated by grain surface reac- 
tions for conditions similar to the knots. The Cloudy models 
confirm this. The grain surface reacti on rate used in Cloudy is 
taken from Cazau x & TielensI (|2002b a). The grain surface for- 
mation rate scales linearly with the integrated projected grain 
surface area. The latter depends critically on the assumed grain 
size distribution as well as the shape of the grains and the grain 
abundance. The integrated grain surface area in our models is 
S = 9.78 X 10"^^ c m^/H. Using other size distributions found 
in the literature (e.g. IWeingartner & Draind 12001.) yielded val- 
ues for S with a range of a factor of 5, with our adopted value 
being towards the upper end of that range. If we assume that the 
grains have a fractal shape rather than spherical, the projected 
grain surface area will be larger than assumed above, yielding 
shorter formation timescales. 

The formation timescales from the Cloudy models can be 
summarized in the fitting formula given in Eq.[T] 



Jh, 



1170 



„0.9 



S-2 



[yr]. 



(1) 



Here «6 is the hydrogen density in the knot in units of 
10* cm"^ and £^21 is the p r ojecte d grain surface area in units of 
10"^' cm^/H. lO'Dell et al] (12007 ) state that the central star has 
exhausted hydrogen-shell nuclear burning one to two thousand 
years ago. Given the observed core densities in the knots, signif- 
icant H2 formation would be possible (though conversion from 
H" to H2 need not be complete) assuming that the knot forma- 
tion was quick after the recombination started and keeping the 



uncertainties stated above in mind. Hence this is currently the 
most plausible explanation for the presence of H2 in the knots. 
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Fig. 2. The H2 mass fraction as a function of the stopping crite- Fig. 5. Same as Fig. [3] but for the standard4 1 
rion. The vertical Hne indicates where T, reaches 4000 K. 





Fig. 3. The standard model of AG04 recreated with Cloudy. 
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Fig. 4. Same as Fig. |2l but for the standard4 model. 

Appendix A: Data reduction 

NGC 6720 was observed by PACS in scan map mode (one scan and one cross 
scan), obtaining maps in the blue (70/im) and red (160/jm) bands. The data were 
proce ssed with the Herschel Interactive data Processing Environment (HIPE, Otl 
I2OIOI) . following the recommended pipeline for these data. This will be explained 
in detail in later papers, but a few deviations and custoinizations are noted here. 
The WCS has an offset of a few arcseconds that was corrected for the compari- 
son of the PACS map to others (using known sources in the field). We used the 
version 3 Hatfield calibration. The signals were converted to Jy using the version 
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Fig. 6. The H2 formation and destruction timescales (solid and 
dotted lines, respectively) for knots with a central density of 10^, 
10^, and 10^ cm"^ just outside the ionized region. 



5 calibration table. We did not remove the cross-talk or coirected for response 
drifts. Cross-talk correction is not part of the pipeline at present. Response drifts 
are unlikely to be a problem for our observations: we find no drift in the cah- 
bration source signals of inore than 1.5% over the duration of the observations. 
Removal of the glitches (cosmic rays) was done in two stages: first from the re- 
gions around our source by filtering along the tiine-dimension, and then from 
the source itself working this tiine on the linage plane. Cleaning the maps of 
the 1/f noise (the name refers to the type of power spectral density the noise 
has) was done using a high-pass filter method. Here a filter passes over the data 
as a function of time, subtracting the median of the data over a specified time 
span (filter width). This allows the high noise frequencies to pass and attenuates 
the lower frequencies. Considering that moving forward in time means moving 
along a spatial direction as the instrument is scanning over the source, it is clear 
that setting a low value for the filter width will rapidly remove varying noise 
but can have the consequence of removing extended emission. We processed the 
data with a vaiiety of filter widths, but for our source it made little dift'erence 
to the photometry or morphology. Also worth noting is that the high-pass filter 
creates artificial negative "sidelobes" around strong sources. To deal with this 
the source was masked out before running the high-pass filter Finally the scan 
and cross scan frames were combined and turned into a map. 

Prior to making the maps presented in this paper the background sources 
were removed with the HIPE sourceExtractorDaophot task. Slight WCS oft'sets 
for the PACS and SPIRE maps were corrected, using background source coor- 
dinates measured in the Calar Alto image. The pixel sizes and beam widths are 
given in Table lB.l1 

We custom-processed the raw Spitzer data obtained from the archive using 
both the GeRT and Mosaicker software to remove some obvious anomalies ap- 
parent in th e pipehne-processed images by following the procedure explained by 
IUetj<2006l) . 
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Table B.l. Aperture photometry of NGC 6720 in various photo- 
metric bands. The apertures are elHptical and the semimajor axis 
is Usted in the second column. The fluxes are listed with their 
measurement and calibration uncertainties, respectively. For the 
Spitzer data no measurement uncertainties were obtained. The 
AKARI data were taken from Izumiura et al. (in prep.). The pixel 
and beam sizes in the Herschel and Spitzer maps are also given. 
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rejected the SWS spectrum as it turned out that the aperture was mainly pointed 
at the central "cavity" and was therefore strongly biased towards the high exci- 
tation region of the nebula. This made aperture correction factors highly depen- 
dent on the ionization stage and hence very uncertain. We did use a re-reduced 
version of the LWS spectrum o f NGC 6720, adopting the aperture correction 
factor listed in iLiu et al] J2004i> . The [Oi] lines were excluded from the mod- 
eling because they were fitted very badly. A plausible explanation is that these 
lines are predominantly formed in the knots. We added two dust continuum flux 
measurements from the LWS spectrum at 43 and 115 fim of 4 X 10""* and 

-5 X 10"" Wcm"^/im"' respectively (aperture correction factors have not yet 
been applied to these values). We also used two radio continuum flux densities 

_at 4850 and 1400 MHz of 3 60 and 440 mj y respectively, which are averages 
of th e data collect ed by Voll mer et alj (20091) . For the ang ular diameter w e used 
76" (Liuetal.'2004). We adopted a distance of 740 pc I O' Defl et alj2009h . The 
model was stopped when the electron temperature dropped below 4000 K. The 
resulting optimized model has the parameters listed in Table [T] Most symbols 
have their usual meaning. T denotes the dust-to-gas mass ratio and e the loga- 
rithmic abundance of an element (e(H) = 12.00). The electron temperature and 
density are averaged over the ionized nebula. We will refer to this model as the 
standard Cloudy model of NGC 6720. 



We did the observations of NGC 6720 in the H2 band by alternating object 
and sky exposures. We reduced the data in the Pixinsight Core package using 
the acquired bias, dark, flatfield and sky frames. Special care was taken to avoid 
artifacts in the object frames caused by the presence of stars in the sky frames. 
A percentile clipping integration was done on the sky frames, grouping them in 
groups of 5 images. By fitting the average signal of the (already bias, dark and 
flatfield corrected) sky frames, a tight rejection of outliers was possible, allow- 
ing the removal of the stars in the images. Then the resulting sky frame was 
subtracted from the object frame acquired in the middle of the five sky frames. 
We corrected the residual background gradients due to sky variation by sub- 
tracting a sky background model built with the DynamicBackgroundExtraction 
module of Pixinsight. F inally this image w as astrometrically calibrated using the 
astrometry.net package jLang et aU2009l) . 



Appendix B: Photometry 

To measure the fluxes from the PACS, SPIRE and Spitzer maps we used an el- 
liptical aperture around NGC 6720. The measured values are given in Table lB.l1 
The PACS data were already in Jy/pix, the SPIRE data were converted from 
Jy/beam to Jy/pix using the conversion ;r(beam/pixel)^/(41n2) where the pixel 
and beam sizes are listed in Table \BA\ The measured fluxes were also multi- 
plied by factors provided by the SPIRE team as the calibration tables did not yet 
include thesfl The Spitzer data were in MJy/sr and were converted to Jy/pix 
via the conversion given in the MIPS instrument handbook (Sect. 4.3). No addi- 
tional corrections were applied (e.g. color corrections). Uncertainties are quoted 
in Table IbD Cahbration uncertainties were taken from the respective instru- 
ment guides or release notes. Measurement errors are difficult to calculate, as for 
these (bolometer) instruments the Poissonian errors are not easily obtained. We 
combined the contribution of the sky noise, the values in the error an'ays which 
the PACS and SPIRE pipelines create, and the variation in map fluxes that dif- 
ferent reasonable pipeline parameter variations gave. For the Spitzer fluxes no 
measurement uncertainties were calculated. 

As the beam size and the pixel scales on the maps are all difl'erent we mea- 
sured all the flux that could be seen from the source down to the sky level, inde- 
pendently for each map. The aperture sizes used are included in Table lB.ll 



Appendix C: The photoionization model 



We used the method described in lvan Hoof & Van de Steena 119991) to create an 
optimized photoionization model of NGC 6720 using a prerelease of version 
CIO. 00 of the photoionization code Cloudy (revision 3862). T he method wa s 
slightly modified in that we used H-Ni model atmosphere s fromlRauchI i2003h . 
As input we used the UV and optical spectrum listed in lLiu et al] 120041) . We 



^ The SPIRE flux calibration multipliers are: 250 /im: 1.02, 350 yum: 
1.05, 500 /im: 0.94. 



